Generalized Pseudopotentials for Higher Partial Wave Scattering 
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We derive a generalized zero-range pseudo-potential applicable to all partial wave solutions to the 
Schrodinger equation based on a delta-shell potential in the limit that the shell radius approaches 
zero. This properly models all higher order multipole moments not accounted for with a monopolar 
delta function at the origin, as used in the familiar Fermi pseudopotential for s-wave scattering. 
By making the strength of the potential energy dependent, we derive self-consistent solutions for 
the entire energy spectrum of the realistic potential. We apply this to study two particles in an 
isotropic harmonic trap, interacting through a central potential, and derive analytic expressions for 
the energy eigenstates and eigenvalues. 

PACS numbers: 34.50.-s 32.80.Pj 34.20.Cf 



The first step in studying the complex physics of a 
many-body system is modeling the fundamental two- 
body interactions. At low energies and for central poten- 
tials, a partial wave expansion reduces the complexity. 
Near zero temperature, s-wave scattering typically dom- 
inates and the true interaction potential can be modeled 
by a contact interaction via the Fermi pseudopotential, 
parameterized by the s-wave scattering length Q| . Such 
a description provides a highly accurate model of the be- 
havior of quantum degenerate gases . For energies out- 
side the Wigner-threshold regime, the Fermi pseudopo- 
tential can be generalized to include an energy-dependent 
scattering length so that the mean-field and trap energy- 
level shifts can be calculated in a self-consistent manner 



Further generalizations are necessary when higher or- 
der partial waves contribute to the interaction. In this 
Letter we derive an energy-dependent pseudopotential 
applicable to I > scattering that captures the criti- 
cal features of both the free and bound-state spectrum 
of the realistic interaction potential. We apply this to 
exactly solve the Schrodinger equation for two particles 
trapped in a harmonic well, interacting through a central 
potential, something typically done only in perturbation 
theory. An example in which higher partial wave scatter- 
ing of trapped particles plays an important role is in the 
physics of degenerate gases of identical fermions where 
the s-wave scattering cross section vanishes |4|. More- 
over, even for bosons, these higher I waves can be reso- 
nantly coupled to the dominant I = scattering due to 
noncentral forces [j| such as the dipolar spin-spin inter- 
action 0, second order spin-orbit interaction Q, and/or 
the forces of an anisotropic trapping potential || 
An example of the latter is the interaction between two 
atoms, each trapped in a separate harmonic well, leading 
to an axially symmetric potential for the relative coor- 
dinate. In previous work, we discovered a new "trap in- 
duced shape resonance" whereby molecular bound states 
are shifted into resonance with trap vibrational levels due 



to the the trap's potential energy [Jjj- Such a resonance 
can have an Z-wave symmetry, and can be modeled by 
our generalized pseudopotential. 

The first attempt to derive a generalized pseudopoten- 
tial was made by K. Huang and C. N. Yang 0, fl2| . 
Given a central force, the true asymptotic wave function 
for each I wave was supposed to follow from a contact 
potential, 
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where ai is the effective scattering length generalized to 
higher partial waves. Recently Roth and Feldmeier (liij 
uncovered difficulties with this pseudopotential noting 
that the mean-field energy shift of interacting fermions 
in a trap is incorrect by a factor of (I + 1)/(2Z + 1). They 
proposed that a distinct effective contact interaction is 
needed to calculate energy level shifts in perturbation 
theory and that the Huang and Yang pseudopotential is 
simply not a proper effective interaction to use in a mean- 
field description of dilute quantum gases Eg. We show 
here that this distinction is unnecessary. Rather, the dis- 
agreement is due to a fundamental problem in Huang's 
original derivation of the pseudopotential E2- Huang 
and Yang incorrectly map the higher order multipolcs 
associated with I > onto a monopolar 6 function at the 
origin. We correct this by employing a <5-shell potential in 
the limit as the shell radius approaches zero. In this limit 
our pseudopotential approaches Huang's original pseu- 
dopotential, but with the correct prefactors thereby giv- 
ing both the correct matrix elements and eigenfunctions. 
Our approach also corrects a different formulation of the 
contact interaction which is used often in the description 
of Rydberg atom collisions E3 • 

In the contact potential construction, one takes the 
asymptotic radial wave function associated with a given 
partial wave Ri(r) = A\ [ji(kr) — tan 5i(k)ni(kr)], valid 
only outside the range of the true potential, and extends 
it to all r. The boundary condition at the origin is set 
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by the zero-range potential, parameterized by the /-wave 
asymptotic phase shift Si(k). Here we assume that the 
realistic potential has a finite range, valid when it falls 
off like 1/r 3 or faster. The scattering phase shift can 
be calculated directly via numerical or analytic solution 
to the Schrodinger equation or may be obtained through 
spectroscopic data. In order to treat the multipole singu- 
larity of the 5 potential at the origin correctly, we write 
the pseudopotential as the limit of a S shell with its radius 
approaching zero, 
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where the operator 0(r) contains the correct prefactors 
and regularization. To derive the correct form of 0(r) 
we solve the radial Schrodinger equation. The familiar 
inside and outside solutions expressed in spherical Bessel 
and Neumann functions are 

R^(r)=B l \ji(kr)} for r < s, (3) 

i?+(r) = Ai [ji(kr) - tan Si(k)ni(kr)] for r > s. (4) 

Requiring continuity of the wave function at r = s hxes 
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where we have used the asymptotic forms of the Bessel 
functions in the limit ks -C 1. Integrating the radial 
equation over the S function gives us a second boundary 
condition. Again taking s -C 1/fc, and using Eq.JSJ), 
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We can fulfill this condition by choosing 
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With the reduced mass fj, and H scaled to one, the pseu- 
dopotential is then 
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Comparing this to the original Huang and Yang pseu- 
dopotential we see that they differ by a factor (i+l)/(2i+ 
1). This occurs because the original derivation ignores 
the inside wave function contribution of weight +1 in 
Eq.0. The (5-shell potential approach circumvents the 
singularity at the origin, allowing one to correctly cap- 
ture higher multipoles. Furthermore, the (5-shell poten- 
tial also enforces the correct ordering of limits, taking 
s — ► as the final step. With this correction, we repro- 
duce the perturbative mean-field energy level shift found 



by Roth and Feldmeier [l3j with a mathematically rigor- 
ous contact potential that also yields the correct asymp- 
totic eigenfunctions. 

In the above form, the (5-shell potential is not Her- 
mitian as the derivative acts solely to the right. The 
regularization we choose is necessary in order to extend 
the domain of the corresponding Hamiltonian to irreg- 
ular functions that diverge as l/r l+1 when r — > 0. Al- 
though this does not cause a problem in most applica- 
tions, in general one must be cautious. In order to make 
the potential Hermitian on the whole domain, including 
both regular and irregular functions, an additional regu- 
larization operator [r l /(2l + 1)!] (d 2l+1 / dr 2l+1 r l+1 ) can 
be added that acts to the left as projector onto the reg- 
ular function subspace. Such dual regularization is cum- 
bersome and so we generally choose to work with only a 
single regularization operator. 

Our form of the (5-shell potential depends on the 
energy-dependent phase shift Si(k) which can usually be 
approximated in the Wigner-threshold regime by a con- 
stant scattering length. We find it more useful here, how- 
ever, to define a fully energy-dependent Z-wave scattering 
length that captures not only corrections due to the ef- 
fective range, but all higher order terms, 
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As in previously studied s-wave case |3j, the general l- 
wave (5-shell potential in our derivation exactly repro- 
duces the correct energy- dependent scattering phase shift 
5i{k) that arises from the true potential and therefore 
exactly reproduces the correct asymptotic wave func- 
tions for all partial waves at all energies. In fact, us- 
ing an energy-dependent scattering length for higher par- 
tial wave scattering has added benefits since the Wigner- 
threshold law may not hold for all I, leading to strong 
energy dependence of the scattering length Eq. (110ft near 
zero energy. For example, for power law potentials of the 
form C n /r n with / > n/2, the phase shift is not propor- 
tional to k 2l+1 but instead behaves as k n [lif. Although 
the generalized scattering length at low energies is not 
constant, the full energy-dependent solution will hold. A 
general breakdown of the pseudopotential approximation 
only occurs in cases where the realistic potential does not 
have a finite range and an outside wave function cannot 
be defined as in Eq.@. 

For two interacting particles in a trap, one solves for 
the discrete eigenvalues of the energy- dependent Hamil- 
tonian derived from an energy-dependent pseudopoten- 
tial using a self-consistent procedure 0. To this end, the 
eigenspectrum of the system is first calculated as a func- 
tion of a constant scattering length, giving E(a). Then 
the effective scattering length is calculated as a function 
of kinetic energy Ek for untrapped scattering states of 
the interaction potential, yielding a(E). Simultaneous 
solutions are then found numerically. This two step pro- 
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cedure allows one to accurately determine the exact shift 
to the energy levels in the trap due to atomic interactions 
when the interaction range and width of the trap wave 
functions are orders of magnitude different, as is typically 
the case. 

For negative energies, k = in is purely imaginary and 
one must analytically continue the scattering length, 



,2H 



tanh [i5i(in)] 



(11) 



Similar to the s-wave case jldj. this analytic continua- 
tion allows us to calculate both the shift in the energy 
spectrum of the trap eigenstates (positive energies) and 
the bound states of the interaction potential (negative 
energies). Consider the radial wave function for negative 
energies, 
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expressed in terms of the spherical Hankel functions of 

(1 2) 

the first and second kind h, ' . Strictly speaking this so- 
lution is only allowed for a normalizable wave functions; 
the true bound states of the 5-shell potential. These oc- 
cur when tanh [iSi(iK)] — 1 since then the coefficient of 

(2) 

the exponentially increasing h\ vanishes. At these en- 
ergies, ai — 1/k. The <5-function bound states are thus 
located at Eg — -(fin) 2 /(2/x) = -h 2 /(2(ia 2 ), just as in 
the s-wave case. The condition for a ^-function bound 
state, tanh [z5;(zk)] = 1, is fulfilled only when the phase 
shift has a pole on the imaginary axis, 8i{in) — —ioo. 
This occurs at each of the negative energies at which the 
S-matrix of the true interaction potential has a pole, i.e., 
at the energies of each of its bound states. The gener- 
alized /-wave pseudopotential with an energy-dependent 
scattering length thus provides an accurate description 
of the entire energy spectrum of the true interaction po- 
tential, bound and scattering, even for /-wave scattering 
lengths with strong energy dependence, i.e. outside the 
Wigner-threshold law regime. 

Our 5-shell approach offers a direct method for obtain- 
ing analytic solutions to a scattering problem by simply 
matching boundary conditions across the S shell as we 
now demonstrate by employing the energy-dependent S 
shell to find all partial wave solutions to the Schrodinger 
equation for two particles in an isotropic harmonic trap 
interacting through a central potential. This is of par- 
ticular interest for application to degenerate quantum 
gases, e. g. two interacting identical fermions. In the 
following, all distances are scaled to characteristic har- 
monic oscillator length zq = H/JjmJ). After separat- 
ing out the center of mass motion, we make the Ansatz 
i? i ± (r) = r l exp(—r 2 /2)w^(r) for the relative coordinate 
radial wave function inside and and outside the shell. The 



radial equation including the scaled trapping potential 
2 then reduces to the Kummer differential equation, 
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zw"(z) + (b — z)w'(z) — aw(z) = 0, in the regions 
where the interaction potential is zero. Independent so- 
lutions of this equation are the confluent hypergeometric 
functions, U(a, b, z) and M(a, b, z), where z — r 2 , a = —u 
and b = 1 + 3/2. The inside solution must be proportional 
to r l exp(—r 2 /2)M(—is,l + 3/2, r 2 ) which behaves regu- 
larly as r l around the origin, whereas the outside solution 
must be proportional to r l exp(—r 2 /2)U{—v,l + 3/2, r 2 ) 
which falls of exponentially for large r. The inside and 
outside solutions are then 

R^{r) = B l r l e- I rM(-is,l + 3/2, r 2 ) for r < s, (13) 
R+(r) = Air l e-4u(-i>,l + 3/2,r 2 ) for r > s. (14) 

We again require continuity of the wave function at r = s. 
In the limit b « 1, 
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where C t = (-l) l T(l + 3/2)T(-v - I - 1/2)/tt and A\ = 
Ai/Ci. Integrating the radial equation over the S function 
gives us again a second boundary condition. Taking s <C 
1/k and using Ea. l(T5|l . the derivatives of the outside and 
inside radial solutions are 
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Applying the operator 0{r) |[SJ| to (|14|) for small s 
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and inserting Eq . Hlfcif) and (|17|l into ||BJ we arrive at the 
implicit eigenvalue equation, 
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This is the general eigenvalue equation for the i-partial 
wave interaction that must be solved self-consistently for 
the energy-dependent ai as described above. For I = 
this reduces to the known s-wave eigenvalue equation 
[r?|. The corresponding wave functions are the inside 
and outside wave functions noted above, where the ratio 
Bi/M is fixed by Ea. (|15[l . For finite shell radius these 
wave functions are in principle numerically normalizable 
unlike solutions obtained with a 5 potential at the origin 
where the unnormalizable solutions diverge as r~(' +1 ) for 
r -> 0. 

In order to verify the accuracy of the higher partial 
wave energy spectrum we choose a spherically symmet- 
ric step potential well with range d and depth Vo as a 
test [see Fig. 1 (a)]. Figures l(b)-l(d) show a comparison 
of the exact eigenspectra and the 5-shell approximation 
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FIG. 1: (a) Step-well test interaction potential with range 
d — 0.4zo and depth Vo = 34.95/kj in a harmonic trap, I 
I eigenstates (dotted lines) and corresponding reduced wave 
function (solid line). The <5-shell solution in the limit of the 
shell radius s — > (dashed line) coincides with the actual 
eigenstates outside the range d. (b)-(d) Comparison between 
exact eigenvalues (solid lines) of the step- well interaction plus 
harmonic trap and pseudopotential eigenvalues (crosses) as a 
function of the range of the well potential and for I = 0, 1, 2 
states. The unshifted interaction bound states (fixed at Et = 
— 2) and trap eigenstates are shown as dashed lines. 



for wells with different finite range d. In particular, we 
chose a well with an Z-wave bound state close to disso- 
ciation to emphasize the accuracy of the approximation 
even in the regime of strongly energy-dependent scat- 
tering lengths ai where the Wigner-threshold law does 
not hold. We find good agreement for relatively large 
ranges d of the well test potential as shown in Fig^ 
The breakdown of the pseudopotential approximation at 
larger ranges is due to the modification of the inter- 
action potential over its finite range by the harmonic 
trap. One can estimate the difference between the en- 
ergy shift with and without this modification for the in- 
teraction bound states in first order perturbation theory, 
AE = (^ s honk 2 /2|V> s hcn) - (ip we n\r 2 /2\ip weU ). Here ip weU 
is the exact bound state associated with the step-well 
potential bound state and ipshell is the bound state wave 
function of the (5-function bound state. For I = 1 and 
I = 2 these two wave functions differ more substantially 
than for s- waves, resulting in a bigger deviation of the 
pseudopotential approximation as the range d becomes 
large [Figs. 1(b)- (d)]. In the case of ultracold collisions 
the energy-dependent pseudopotential will therefore be a 
good approximation as long as the characteristic interac- 
tion length scale of the van der Waals interaction is much 
smaller that the characteristic length scale of the trap zq 



In summary, we have derived a generalized zero-range 
pseudopotential for higher partial wave interactions that 
captures both the scattering solutions and bound-state 
spectrum self-consistently. By employing a limiting pro- 
cedure on a finite radius (5-shell potential, we provided a 
rigorous correction to the long standing error in Huang's 
and Yang's pseudopotential. The pseudopotential offers 
a direct method to analytically solve the Schrodinger 
equation, as demonstrated for the case of interacting 
trapped atoms, where we derived the higher partial wave 
energy spectrum and obtained normalizable eigenfunc- 
tions. This is of special interest for degenerate gases 
of identical fermions where I — 1 scattering is the pri- 
mary contribution to the interaction and also for Bose 
systems where noncentral forces play an important role. 
Our accurate modeling of the interaction and the analyt- 
ical calculation of the eigenenergies should provide new 
avenues for studying degenerate gases of interacting ul- 
tracold atoms in tightly confining traps 0, such as in 
optical lattices |18|. Beyond its application to many- 
body problems the (5-shell pseudopotential is also useful 
for modeling controlled collisions, which could play an 
important role in quantum information processing. 

After completing this paper, we learned of recent theo- 
retical work by Kanjilal and Blume |9( in which the I — 1 
special case of Eq. (|18fl has been derived and applied to 
ID and 3D confined fermions, and of work by Peach et al. 
|l9l | in which Eq. I |18| l has been derived using a quantum- 
defect theory approach. 

We thank Paul Julienne, Carl Williams, and 
Eite Tiesinga for very helpful discussions. This 
work was partly supported by NSA/ARDA Contract 
No. DAAD19-01-1-0648 and by the ONR Contract 
No. N00014-03-1-0508. 
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